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Abstract 

We present studies of the photoexcited quasiparticle dynamics in MgE>2 where, using femtosecond 
optical techniques, Cooper pair breaking dynamics (PBD) have been temporally resolved for the 
first time. The PBD are strongly temperature and photoexcitation intensity dependent. Analysis 
of the PBD using the Rothwarf- Taylor equations suggests that the anomalous PBD arises from 
the fact that in MgB2 photoexcitation is initially followed by energy relaxation to high frequency 
phonons instead of, as commonly assumed, e-e thermalization. Furthermore, the bare quasiparticle 
recombination rate and the probability for pair-breaking by phonons have been determined. 
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In recent years, femtosecond real-time spectroscopy has been shown to present an excel- 
lent experimental alternative for studying temperature (T) dependent changes in the low 
energy electronic structure of strongly correlated electron systems Q Q li Q] In 
these experiments, a femtosecond laser pump pulse excites electron-hole pairs via an inter- 
band transition in the material. In a process which is similar in many materials including 
metals and superconductors (SC), these hot carriers rapidly thermalize via electron-electron 
(e — e) and electron-phonon (e —ph) collisions reaching states near the Fermi energy within 
10-100 fs. The subsequent relaxation and recombination dynamics (strongly affected by 
the opening of the superconducting or charge density wavep] gap), are ob- 

served either by measuring photoinduced (PI) changes in reflectivity or transmission at 

|, or by measuring conductivity dynamics at terahertz (THz) 
While higli^emrjerature superconductors have been extensively studied over 
1 



optical frequencies [1, 
frequencies 



the last decade or so 
is fairly limited P-llfl 



^| , the data on more conventional BCS type superconductors 
. To our knowledge, no systematic study of photoexcited quasiparticle 
(QP) dynamics with femtosecond resolution has been performed on low-T c superconductors. 

In this Letter, we present the first femtosecond time-resolved study of photoexcited carrier 
dynamics in the recently discovered superconductor MgB 2 , utilizing both optical pump — 
THz probe (OPTP) and optical pump — optical probe (OPOP) techniques. The aim of this 
study is to elucidate the photoexcited carrier dynamics in MgB 2 and in SC's in general, to 
compare the dynamics to those in cuprates, and to obtain new and complementary knowledge 
of the electronic structure of MgB 2 . The discovery of bulk superconductivity below 39 
K in MgBojlll| has generated a great deal of excitement since T c is higher by nearly a 
factor of two in comparison to other previously known simple intermetallic superconductors. 
While the observation of a significant boron isotope effect 3], and the spin-singlet nature 
of the pairing 
experimental 
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13 1 indicate conventional phonon mediated BCS pairing in MgB 2 , several 

™ff theore ' ical studies k™ ^ — of 

two distinct energy gapsjljj. 

Since the SC gaps in MgB 2 lie in the THz range^Jj, OPTP spectroscopy is well-suited 
to study the dynamics of PI quasiparticles. In particular, the Cooper pair-breaking dynam- 
ics (PBD) have been time-resolved for the first time. The PBD are strongly temperature 
and photoexcitation fluence dependent. This is attributed to, following photoexcitation, an 
initially strong relaxation to high frequency phonons, and supported by detailed analysis in 
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21| . The SC state recovery dynamics, 



terms of a phenomenological Rot hwarf- Taylor model 
on the other hand, proceed on the timescale of hundreds of picoseconds. Systematic stud- 
ies as a function of temperature, excitation intensity, and film thickness suggest that pair 
recovery is governed by the anharmonic decay of 2A acoustic phonons. 

In time-domain THz spectroscopy a nearly single-cycle electric field transient, containing 
Fourier components from ~ 100 GHz to several THz is generated via optical rectification of 
150 fs optical pulses in a ZnTe crystal. The THz electric field transmitted through a sample, 



E sam (t), is detected using the Pockels effect in ZnTe|5j,|7[. A measurement of the transmitted 
electric field is also made using a suitable reference, E re f(t), which, for these experiments, 
is a blank sapphire substrate. Dividing the Fourier transforms of the time domain sample 
and reference data gives the complex transmissivity T(uS) = E sam (uS) / E re f(uj) . The real and 
imaginary conductivities (<J r (u) , <Ji(u))) of the film are determined using the appropriate 
complex Fresnel equation, without the need for Kramers- Kronig analysis [22]. For the OPTP 
experiments, the induced change in the transmitted electric field, AE sam (t), is measured as 
a function of time-delay with respect to an optical excitation pulse. From AE s<m {t}, it is 



,0. 



then possible to determine the PI change of cr(a>) with picosecond (ps) resolution [71. \2c 

The OPTP experiments were performed on 80 nm and 100 nm MgB 2 thin films (T c ~ 34 
K) on sapphire 251 ]. whereas the films used in OPOP experiments had thicknesses of 300 and 
400 nm, and T c = 39 K (AT C ~ 0.15K). Further details of the fi m growth 3l and a detailed 
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The photoexcitation 



description of experimental techniques are given elsewhere 

fluence F ranged from 0.1 — 5 /iJ/cm 2 , corresponding 2f| to an absorbed energy density 
Q = 2 — 110 /zeV/unit cell. For comparison, Q corresponding to complete destruction of SC 
state was found to be ~ 110 //eV/unit cell, in agreement with the condensation energy in 
the BCS limit £ | . 

Figure 1 shows <Ji{oj) and a r (u) at several time delays aft er p hotoexcitation with a 150 
fs pulse with F ~ 3 /iJ/cm 2 . Since Oi provides a direct probe [22] of the condensate density, 
n s , by measuring PI changes in Oi direct information of n s dynamics can be extracted. On 
the other hand, the increase and subsequent recovery of a r (u) corresponds to an initial 
increase in the number of QPs followed by their recombination |5j. The decrease of (Ti(u) 
— corresponding to a reduction of n s — occurs on the timescale of several ps, followed by 
recombination dynamics on the timescale of several hundred ps. To obtain a more detailed 
time evolution of the PI changes in a(u), we utilize the fact that the induced changes in 
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the electric field transient, AE sam , are mainly due to a phase shift of the electric field - see 
inset to Fig. la). The origin of the phase shift is the so-called kinetic inductance due to SC 
pairing, i.e. the conductivity due to the SC pairs is purely imaginary with al/w dependence 
resulting in the overall phase shift in E sam (t) below T c . When F is smaller than the fluence 
corresponding to destruction of the SC state, AE sam (t = t Q ) (t is a fixed point of E sam (t) 
which, for these experiments, was at the point of maximum time derivative of the electric 
field - indicated by the arrow in inset to Fig.l) is proportional to the PI conductivity Act 
(both AcTj and Aa r have the same dynamics - see inset to Fig. lb). Therefore by measuring 
AE sam (t = to) while scanning the pump line, the PI conductivity dynamics Aa(t) can be 
extracted rapidly (avoiding the signal drifts associated with the long term laser stability) 
and with much higher temporal resolution. 

Fig. 2a) presents the induced conductivity dynamics as a function of T (F « 1 /ij/cm 2 ). 
The relaxation time tr, obtained by fitting the data to exp(—t/Tji), shows a pronounced T- 
dependence, plotted by circles in the inset to Fig. 2a). Upon increasing T, tr first decreases, 
reaches a minimum, followed by a quasi divergence as T c is approached. 

Similar results are obtained from the OPOP experiments - Fig. 2b), where PI changes 
in reflectivity (AR/R) at optical frequencies are measured (20 fs pulses at Tuo ~ 1.54 eV 
have been used as a source of both pump and probe pulses). The normal state dynamics 2]| 
are characterized by a resolution limited risetime (40 fs), followed by a two-exponential 
decay, with decay times T\ ~ 0.15 ps and r 2 ~ 3.5 ps (Fig. 2b)). After the initial ps 
dynamics, AR/R changes sign and the recovery dynamics proceed on a timescale longer 
than 1 ns. Below T c , however, an additional sub- nanosecond response becomes evident, 
whose amplitude and recovery dynamics is strongly T-dependent }28| . In particular, the T- 
dependence of tr, plotted in the inset to Fig. 2b), displays the same T-dependence as the 
OPTP data. 

The SC state recovery dynamics do not show any dependence on F in the range 
0.1 < F < 5 /iJ/cm 2 . This suggests that the SC recovery dynamics are governed by the 
phonon-bottleneck mechanism, proposed by Rothwarf and Taylor [2jj (biparticle recombina- 
tion should be intensity dependent). In this case recovery dynamics a re g overned by the 
lifetime of u > 2 A phonons, which are in thermal equilibrium with QPs. |2lJ Since in MgB 2 
2A is much smaller than the energy of the lowest optical phonon mode (~ 40 meV) Q, the 
SC recovery is governed by the decay of acoustic phonons. The decay of u > 2A phonon 
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population is governed either by anharmonic decay to u < 2 A phonons or escape of u > 2A 



phonons to the substrate 



281 ] . The fact that no change in tr is observed when compar- 



ing the data taken on 80, 100 and 400 nm films suggests that in MgB 2 the predominant 
mechanism that governs the SC condensate recovery is anharmonic phonon decay. This as- 
signment is further supported by the observed T-dependence of tr, which near T c ^| shows 
tr oc 1/A(T), and the estimated anharmonic decay timej^l of the 10 meV longitudinal 
acoustic phonon in MgB2 of ~ 1 ns 26]. 

When discussing the OPOP data we should mention the peculiar ps normal state dy- 



namics, present also in the SC state. In meta 
terms of the two-temperature model (TTM) 



Sj ps QP dynamics are usually interpreted in 
3lL l32j ] . Here the assumption that the e-e 



scattering is much faster than the e-ph scattering leads to the description of the PI tran- 
sient in terms of the time evolution of the electronic temperature T e , i.e. AR(t) = dR/dT 
AT e (t). The dynamics are due to e-ph thermalization which is proportional to the e-ph 
coupling constant A^l^. For MgB2, however, the ps dynamics are inconsistent with the 
TTM since: i) the AR/R transient changes sign with timej^J, ii) dR/dT at 1.5 eV is 
negative pS], so AR/R < in the TTM, contrary to what is observed, iii) the expected 
e-ph Nation WQ is r ep ?s 20 fs, which is much shorter than the experimental 
values of either T\ or T2, and iv) the initial dynamics do not change upon cooling below T c , 
where the bottleneck in the relaxation due to the presence of the SC gap should strongly 
affect the recovery dynamics [ll. I^. The above arguments suggest that the origin of the ps 
dynamics measured in OPOP experiments is not associated with e-ph thermalization, and 
indicates that the assumption of e-e thermalization being much faster than the e-ph scatter- 
ing is invalid in MgB 2 . The coupling between electrons and high frequency optical phonons 
in MgB2 is very strong ( in p articular the coupling of o band electrons to the E2 9 phonon 
mode at 60-80 meV) [3, 13 • Thus, it is quite possible that the situation in MgB 2 is 
reversed (or at least that the two timescales are comparable). In this scenario the initial 
relaxation of photoexcited electrons proceeds via emission of high frequency (60-80 meV) 
optical phonons which only subsequently release their energy to the electron system via 
phonon-electron scattering, and to low energy phonons via anharmonic decay. Therefore we 
argue that the ps dynamics observed in OPOP experiments is due to the energy relaxation 
of the optical phonon population, rather than e-ph thermalization. 

Finally, let us discuss the rise-time dynamics, reflecting the Cooper-pair breaking pro- 
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cesses (i.e. the initial reduction of n s ). Fig. 2a) and the inset to Fig. lb) clearly show a 
finite risetime in the induced change in conductivity, indicating that it takes some time for 
the completion of pair-breaking following optical excitation. Furthermore, Fig. 2a) shows 
that the PBD are T-dependent (the PBD becomes faster as T is increased). Also, the PBD 
depend on the photoexcitation intensity. Fig. 3a) shows the early time Aa(t) taken at 7K 
for different fluences. While the solid symbols represent the data obtained by measuring 
AE sam (t), the open symbols represent the induced change in the conductivity obtained di- 
rectly using the two-dimensional scanning technique. The agreement between the two data 
sets clearly shows that the phase changes accurately reveal the condensate dynamics. This 
is particularly important, since measuring AE sam (t) enables the study of condensate dy- 
namics with sub-picosecond resolution (limited by optical pulsewidths to ~ 0.3 ps), while 
the two-dimensional scanning technique is limited by the THz pulse width (~ 2 ps) |24j |. 

In view of the normal state OPOP data, which suggest that high energy electrons initially 
release their energy (or at least a significant portion of it) via emission of high frequency 
phonons (which only subsequently break Cooper pairs), there is a natural explanation of 
the observed fluence and T-dependence of the PBD in MgB2 j28j] . To show this, we consider 
the phenomenological Rothwarf- Taylor model[21], where the dynamics of the QP and high 
frequency (u > 2 A) phonon densities, n and N, are described by a set of two coupled dif- 
ferential equations |21 1 . Since the SC recovery dynamics proceed on a much longer timescale 
than the PBD, the term describing the loss of u > 2A phonons by processes other than pair 
excitation can be neglected, and the equations are given by: 

dn/dt = pN - Rn 2 ; dN/dt = ^[Rn 2 - /3N] (1) 

Here R is the bare quasiparticle recombination rate and f3 is the probability for pair-breaking 
by phonons 21}. With the initial condition that after photoexcitation (and initial sub- 
picosecond e-e and e-ph dynamics) the QP and high frequency phonon densities are n 
and N , the subsequent time evolution of n is given by 
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28] 

(2) 



4 2r rl- Kexp(-tf3/r 
where K and r are dimensionless parameters determined by the initial conditions: 



Eq. [2] has three distinct regimes, depicted in inset to Fig. 3a). K = corresponds to the 
stationary solution, when no and No after photoexcitation are already in quasi-equilibrium 
(at a somewhat elevated T), i.e. Rhq = (3N Q , and n(t) is a step function. The regime 
< K < 1 corresponds to the situation when, following excitation, the number of QPs is 
higher than the quasi-equilibrium value, while —1 < K < represents the opposite situation 
when the initial T of u > 2A phonons is higher than the quasi-equilibrium one. As Fig. 
3 shows, the latter situation seems to be realized in the case of MgB 2 , consistent with the 
analysis of the OPOP data. 

At low T, when the density of thermally excited QPs and uj > 2 A phonons is exponentially 
small, the change in conductivity is proportional to the PI quasiparticle density. Therefore 
we can the conductivity dynamics using Eq. ({%))■ Best fits to the data taken at various 

intensities are plotted by solid lines in Fig. 3a) showing extremely good agreement with the 
data. The two fitting parameters t//3 and K are plotted in panels b) and c). 

Importantly, the microscopic parameters R and (3 can be extracted by fitting the t/(3 
and K vs. the absorbed energy density Q to Eqs.(jBJ). We assume that photoexcitation 
with Q creates no = pQ/A QPs (with energy A) and Nq = (1 — p)Q/2A high frequency 
2A phonons, where p is the portion of Q that initially goes to QPs. A best fit to the two 
data sets (dashed lines) is obtained when p 6% (further supporting OPOP data analysis) 
giving /T 1 = 15 ± 2 ps, R = 100 ± 30 ps^unit cell" 1 . 

In conclusion, we have presented the first femtosecond studies of photoexcited carrier 
dynamics in MgB2 using both OPTP and OPOP techniques. The T-dependence of the 
SC state recovery dynamics have a similar behavior as for the cuprates^, Q, 0]. The main 
difference is that in the cuprates tr is a factor of 100 shorter, which we attribute to their 
larger gap value, causing the dynamics (and tr) to be governed by the lifetime of optical 
phonons, instead of acoustic phonons as in MgB 2 j3|. Furthermore, we have presented 
the first observation of time-resolved PBD. The PBD are found to be temperature and 
fluence dependent, which is attributed to the initial creation of high frequency phonons, 
rather than QPs, and supported by detailed analysis of the PBD in terms of the Rothwarf- 
Taylor equationsj2l|. Importantly, similar physics may be responsible for the longer risetimes 
observed in single layer cuprates (~ 500 fs) |4( in comparison to double layered ones (~ 0.1 

ps)Q,flS. 
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Figure Captions: 

Figure 1: The a) imaginary and b) real conductivity as a function of frequency shown 
at various time delays following excitation with a fluence F ^ 3 /zJ/cm 2 at 7K. Inset to a): 
E sa m{t) & t 7K and 35K. Inset to b): the time evolution of a r and <7j taken at v — 0.8 THz. 

Figure 2: a) The PI change in conductivity (Aa oc AE sam ) as a function of temperature, 
b) The PI reflectivity traces at various temperatures plotted on a semi-log scale. The normal 
state trace (90K) is vertically shifted for clarity. The T-dependence of the superconducting 
state recovery time, tr, determined by single exponential fits to the data (solid lines) is 
presented in the insets. 

Figure 3: a) The rise time dynamics at 7 K taken at various F in /iJ/cm 2 (for presenta- 
tion purpose all traces have been shifted vertically by 0.1). Solid circles represent the data 
obtained by scanning the pump line while measuring AE(t = t ), while the open circles 
correspond to Aa measured directly. Solid lines are fits to the data using Eq.(2), with the 
intensity F (absorbed energy fi) dependence of r and K plotted in panels b) and c). Dashed 
lines in b) and c) represent the best fit to Eqs.(3). Inset to panel a): Solutions of Eq.(2) for 
different K. 
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